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Abstract. We consider the zero sets of systeras of m random polynomials of degree N 
in m complex variables, and we give asymptotic formulas for the random variables given by 
summing a smooth test function over Z^. Our asymptotic formulas show that the variances 
for these smooth statistics have the growth iV™~^. We also prove analogues for the integrals 
of smooth test forms over the subvarieties defined by fc < m random polynomials. Such 
linear statistics of random zero sets are smooth analogues of the random variables given by 
counting the number of zeros in an open set, which we proved elsewhere to have variances of 
order N™~^^^. We use the variance asymptotics and off-diagonal estimates of Szego kernels 
to extend an asymptotic normality result of Sodin-Tsirelson to the case of smooth linear 
statistics for zero sets of codimension one in any dimension m. 



1. Introduction 

This article is concerned with zero sets of systems of Gaussian random polynomials (or 
more generally, of sections of a positive holomorphic line bundle over a compact Kahler 
manifold Mm) as the degree N — >■ oo. One of the most fundamental statistical quantities is 
the number Mj!}{pi , . . . , Pm) of zeros in a bounded open set U C of a system {pi, . . . , p^} 
of m independent Gaussian random polynomials. The expected value of this random variable 
was shown in [SZlj to be the integral of the Kahler volume form over U (times a universal 
constant). In a recent article [SZ4], we gave an asymptotic formula for the variance of this 
random variable. We also give analogous results for the volume of the simultaneous zero 
set of < m polynomials or sections. In this article we apply the methods of |SZ4] to the 
analogous 'smooth linear statistics', i.e. the sum (or integral) of a smooth test function 
over the zeros of a system of random polynomials. Such smooth linear statistics arise as 
smooth approximations for discontinuous random variables such as and also arise in 
a number of other problems (see the discussion in Our main results give asymptotic 
formulas for the variance of these smooth linear statistics. (Mean value asymptotics for 
these random variables were given in |SZ1] .) As may be expected, the variances are of much 
lower order in the degree N than in the non-smooth case. Further, we prove that in the 
codimension one case, the smooth linear statistics are asymptotically normal, extending a 
result of Sodin-Tsirelson |STj . 

To state our results precisely, we need some notation and background. We let (L, h) ^ M 
be a positively curved Hermitian holomorphic line bundle over a compact complex manifold 
M, and we give H^{M, L^) the Hermitian Gaussian measure induced by h and the Kahler 
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form ijj = |9/i (see Definition 12 .ip . For m independent random sections Sj^ G H^{M, L^), 
1 < j < m, the number of simultaneous zeros of the sections in a smooth domain U G M is 
given by 

^^{s^, ...,sZ):=nzeU:sn^) = --- = s^^z) = 0}. 

In [ SZlj . we proved that the expected value E(A/'^) of the random variable A/"^ has the 
asymptotics 

N" 



E«) 



vr'' 



uj"' + 0{N"'-') 



u 



and in |SZ4j . the variance of the random variable is shown to have the asymptotics, 
Var(Ar^) =iV™-V2 L VoW_i(af/)+0(iV-i+^ 



(2) 



where z/^m is a universal positive constant ( |SZ4t Theorem 1.1]). Analogous results are proved 
for the volumes of zero sets of A; < m — 1 independent random sections G H^{M,L'^), 
^ ^ j ^ k: in this case, 



VaT{Vohm-2k[z,N_^N n u]) = iv^'^-^-i/s L , voi2™-i(9t/) + o(Ar-5+-) 



, (3) 



where i^mk is a universal constant ( |SZ4l Theorem 1.4]); in particular, z/^i = 
More generally, the domain could be piecewise smooth without cusps. 

In this article, we are interested in the smooth analogue of A/"^ where we integrate a 
smooth test function rather than the characteristic function of a smooth domain over the 
zero set. Given a test function (f G V{M), we consider the random variable 



s^(z)=-=s^(z)=0 



When the system is not full, we define 



m—k,m—k 



The expected value of (Z^jv^ gjv,*/?) is given by (see f|T4l) ) 



(M) 



J M 



(4) 



The main result of this article is an asymptotic formula for the variance: 



Theorem 1.1. Let {L,h) be a positive Hermitian holomorphic line bundle over a compact 
Kdhler manifold {M,uj), where uo = |6/i, and let 1 < k < m. We give H^{M, L^) the 
Hermitian Gaussian measure induced by h,uj (see Definition \2.1\) . 

Let ip be a real {m — k,m — k)-form on M with coefficients. Then for independent 
random sections , ■ ■ ■ G H^{M, L^), we have 



2k-m-2 



M 
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where ^Im is the volume form on M , and Bmk is a universal Hermitian form on the bundle 
T*m-k+i,m-k+i^^y ly/ien k = l, we have Bmi (f^M, f^n) = \ f\\ and hence 

Var(Z,^, ^) = iV— Wdd^Wl, + 0{N-^ 

In particular, for the complex curve case m = 1, we note that \ddip\ = ^\Alp\, and thus 

■C(3) 



Var(Z,iv,(/?) = N-^ 



IQtt 



(5) 



The leading term in ([5]) was obtained by Sodin and Tsirelson |ST] for the case of random 
polynomials G if°(CP^, 0{N)) as well as for random holomorphic functions on C and on 
the disk. (The constant was given in a private communication from M. Sodin.) 

Here we say that B^k is universal if there exists a Hermitian inner product B^i^ on 
rp*m-k+i,m-k+i^^rn^^ independent of M and L, such that for all w G M and all unitary 
transformations r : Tq{C"^) T^{M), we have B^kiw) = t^B^^^. The global inner prod- 
uct {^.Tp) = Jj^,j Bmk {ddip,ddip) Qm is certainly positive semi-definite on p™-'^'™-'^(M), 
since the variance is nonnegative. We believe that, in fact, B'^j^ is positive definite on 
rp*m- +i,m- +^^(^my This follows for k = 1 from the above formula for Bmi', one should be 
able to verify positivity for /c > 1 by using the expansion (!27|) in the proof of Theorem 11.11 
to compute a precise formula for -B^^. 

Thus the variance of the 'smooth statistic' (Z^iv^ ,,jv, (f) is of lower order than the variance 
of the number and volume statistics given by ([2])-([3]), as expected. In view of ([1]), it is also 
self-averaging in the sense that its fluctuations are of smaller order than its typical values. 

An application of our methods is an extension of the Sodin-Tsirelson |ST] asymptotic nor- 
mality result for smooth statistics to general one-dimensional ensembles and to codimension 
one zero sets in higher dimensions: 

Theorem 1.2. Let (L, h) — > (M, u) be as in Theorem{Tj\ and give H^{M, L^) the Hermitian 
Gaussian measure induced by h,uj. Let ip be a real (m — l,m — l)-form on M with coef- 
ficients, such that ddif ^ 0. Then for random sections in H^{M,L'^), the distributions 
of the random variables 

converge weakly to the standard Gaussian distribution Af{0, 1) as N oo. 

Sodin and Tsirelson |ST] obtained the asymptotics of Theorem 1 1.21 for random functions on 
C, CP^, and the disk. The proof of Theorem 1 1.2 1 is a relatively straightforward application of 
the fundamental Szego kernel asymptotics underlying Theorem 11.11 to the argument in [STj. 
(One easily sees that the random variable {Z^n, (p) is constant for all if ddip = 0.) 

Substituting the values of the expectation and variance of {Z^n, ip) from (jlj) and Theorem 
II. H respectively, we have: 

Corollary 1.3. With the same notation and hypotheses as in Theorem \1.2i the distribu- 
tions of the random variables N'^/'^{ZgN — ^uj,(p) converge weakly to A/'(0, a/k^ ||99(y9||2) as 

N — i> oo, where = 
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Here, A/'(0, cr) denotes the (real) Gaussian distribution of mean zero and variance o"^. 

We now summarize the key ideas in the proofs in |SZ4j and in this paper. The variance in 
Theorem ll.il as well as the number and volume variances in |SZ4] . can be expressed in terms 
of the variance currents Var(Z3jv^ giv) of the random currents Z^iv ^jv. In joint work with 
P. Bleher in 2000 [BSZlj . we introduced a bipotential Qn for the 'pair correlation function' 
K21 of the volume density of zeros of random sections in H^{M, L^)\ this bipotential satisfies 

The bipotential Qat is a universal function of the normalized Szego kernel (see (fT6l) and (!22|) ). 
Sodin and Tsirelson |ST] obtained a variance formula as well as asymptotic normality for 
zeros of certain model one-dimensional random holomorphic functions by implicitly using 
this bipotential. 

In |SZ4j . we showed that Qat is actually a 'pluri-bipotentiaF for the co dimension-one 
variance current; i.e., 

iidd)^ {idd)yj Qn{z, w) = Var(Z^iv) . (6) 

We further found a formula (Theorem 13. ip for the higher codimension variance current 
Var(Z3jv^ ,,jv) in terms of Qtv and its derivatives of order < 4. We then applied the off- 
diagonal asymptotics of the Szego kernel IIn{z,w) in |SZ2] to obtain asymptotics of the 
bipotential Qn{z, w) and then of the number variance ([2]) as well as the volume variance ([3]). 

In this paper, we begin by reviewing basic facts about the Szego kernel and summarizing 
the asymptotics from |SZ4j of the bipotential Qat as ^ 00 as well as d{z,w) 0. To 
illustrate our ideas, we apply these asymptotics to to derive the codimension one formula 
(i.e., the case = 1) of Theorem ll.il We then prove in §3.21 a slight modification (Corollary 
13. 3p of the formula for the higher codimension variance, which we use in §3.31 to prove 
Theorem II. 1[ In §1] we apply our Szego kernel asymptotics to prove Theorem 11.21 Finally, 
we formulate some related open problems in ^ 

2. Background 
In this section we summarize results from [SZ4j used in this paper. 

We let (L, h) be a Hermitian holomorphic line bundle over a compact Kahler manifold M. 
We consider a local holomorphic frame over a trivializing chart U. If s = fei is a section 
of L over U, its Hermitian norm is given by ||s(z)||/i = a{z)~2\f[z)\ where 

aiz) = \\eLiz)\\f . (7) 

The curvature form of [L, h) is given locally by 

Qh = ddloga , 

and the Chern form Ci(L, h) is given by 

c,{L,h) = = log a . (8) 

The current of integration Zs over the zeros of a section s G H^{M, L) is then given by the 
Poincare-Lelong formula, 
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where the second equahty is a consequence of ©-([HI). 

We now assume that the Hermitian metric h has strictly positive curvature and we give 
M the Kahler form 

uj='-Qh = 'i^cimh) . (10) 

Next we describe the natural Gaussian probability measures on the spaces H^{M, L^) of 
holomorphic sections of tensor powers = L®^ of the line bundle L: 

Definition 2.1. Let {L,h) —>■ {M,uj) be as above, and let denote the Hermitian metric 
on induced by h. We give H^{M,L^) the inner product induced by the Kahler form uj 
and the Hermitian metric : 



Jm ml 



(11) 



The Hermitian Gaussian measure on H^{M, L^) is the complex Gaussian probability mea- 
sure 7Ar induced by the inner product ( ITTl) .- 

dlN{.s) = — e~l^l'dc , = X] ^J'^f ' 

i=i 

where {Sf, . . . , S^^} is an orthonormal basis for H^{M, L^). It is of course independent of 
the choice of orthonormal basis. 

The Gaussian ensembles (i/°(M, L^), 7^) were also studied in jSZTl [SZ2l[BSZn [BSZ2] : 
for the case of polynomials in one variable, they become the SU(2) ensembles in [BBLt IHat 
INVt IZhj : for polynomials in m complex variables, they are the SU(m + 1) ensembles (see, 

e.g., [SzHEi]). 

We consider the diagonal Szego kernels 

n^(z,z) :=$^||5f(^)||^.(,), 
i=i 

where the are as in the above definition. It follows from the leading terms of the asymp- 
totic expansion of the diagonal Szego kernel of |Cal [Ti| [Ze] that 

IIn{z,z) = \n^{1 + 0{N~')). (12) 

The expected value of the zero current of a random holomorphic section in H'^[M, L^) is 
given by the basic formula 

i - N 
EZsN = —dd\ogUN{z,z) + —iu , (13) 

and the expected values of simultaneous zero currents are given by 

E{Z,._,.) = [E(Z,.)]"' = Q-dd\ogUr,{z,z) + ^uy = ^ + 0{N'~') , (14) 

for 1 < /c < m (see |SZ1[ ISZ4] ). The final equality of f|T4l) is a consequence of the asymptotic 
formula ([12]). 
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For the variance asymptotics, we need the properties of the off-diagonal Szego kerneh 



IIn{z,w)\ :-- 



dN 



w 



(15) 



In particular, our variance formulas are expressed in terms of the normalized Szego kernel 



Pn{z,w) :-- 



\Un(z,w) 



n 



N 



'z,z)2IIn{w,w)2 



(16) 



In |SZ4] . we used the off-diagonal asymptotics for IItv from |SZ2] to provide the off-diagonal 
estimates given below for the normalized Szego kernel Pn- As in |BSZlt ISZlt ISZ2t ISZ4] ). we 
obtained these asymptotics by identifying the line bundle Szego kernel |n7v(^, w) \ of (fT5i) with 
the absolute value of a scalar Szego kernel IIn{x, y) on the unit circle bundle X C L^^ — > M 
associated to the Hermitian metric h. 

Our estimates are of two types: (1) 'near-diagonal' asymptotics (Proposition 12.31) for 
P]\f{z,w) where the distance d{z,w) between z and w satisfies an upper bound d{z,w) < 

1 /2 

6 (-2^) (b G M+); (2) 'far-off-diagonal' asymptotics (Proposition 12.21) where d{z,w) > 



Proposition 2.2. [gZ4l Prop. 2.6] Let {L,h) {M,uj) be as in TheoremUJl and let 

Pn{z,w) be the normalized Szego kernel for H^{M,L^) given by (1161) . For b > y/j + 2k, 
j, k > 0, we have 



Pn{z,w) = 0{N ^) uniformly for d{z,w)>b 



logN 



N 



Here, V-' stands for the j-th covariant derivative. The normalized Szego kernel P/v also 
satisfies Gaussian decay estimates valid very close to the diagonal. To give this estimate, we 
consider a local normal coordinate chart p : U,Zo ^ C™', centered at a point zq & U C M, 
and we write, by abuse of notation, 

Pn{zo + u,Zo + v) := Pn{p'^{u),p'^{v)) . 

Proposition 2.3. |SZ4[ Prop. 2.7-2.8] Let Pn{z,w) be as in Proposition \2/4 and let zq e 
M. For e,b > 0, there are constants Cj = Cj{M,e,b), j > 2, independent of the point zq, 
such that 

P^ (zo + ^,zo + ^)= e--2\^-^\'[l + RNi%v)] , 

where 

\RNiu,v)\ < ^\u-v\'^N~^/^+' , \VRn{u)\ <C2\u-v\N-^/^+', 
\V^Rn{u, v)\ < Cj j>2, 

for \u\ + \v\ < fei/IogTV. 
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2.1. The pluri-bipotential for the variance. For random codimension k zeros, we have 
the variance current of Z„n „jv: 

Var(Z,^,...,,^) = E{Z,._,.MZ,M_,.)-E{Z^M_^.)ME{Z,._^.) e V"''^'\M x M). 

(17) 

The variance for the 'smooth zero statistics' is given by: 

YaT{Z^M_,N,^) = (YRr{Z,N_,N), <^ K v^) . (18) 

Here we write 

RMS = ttIRAtt^S e V'P+i{M X M) , for i? G V'p{M), S G V"^{M) , 

where tti, : M x M — ^ M are the projections to the first and second factors, respectively. 
For a current T on M x M , we shall write 

dT = diT + d2T, di = ydz, ^ , d2 = ydwj^, 

OZj ^-^ OWj 

where zi, . . . , 2;^ are local coordinates on the first factor, and wi, . . . , Wm are local coordinates 
on the second factor of MxM. In particular, di{RmS) = {dR)MS and d2{RMS) = Rm{dS). 
We similarly write 

BT = diT + d2T . 

In |SZ4] . we constructed a pluri-bipotential for the variance current in codimension one, 
i.e. a function Qn & L^{M x M) such that 

Var(Zs]v) = -didid2d2QN = {idd)^ (idd)^ Qn{z,w) . (19) 

To describe our pluri-bipotential Qn{z,w), we define the function 



n=l 



Alternatively, 



G{e-') = -— / log(l - e-'n ds, A > . (21) 



Theorem 2.4. [SZil Theorem 3.1] Let {L,h) {M,uj) be as m TheoremUJi Let Qn 
M X M — > [0, +oo) be the function given by 

l_ fP^i'^^)' log(l - s 



Qn{z, w) = G{Pn{z, w)) = -— I ds , (22) 



where Pn{z,w) is the normalized Szego kernel given by ( fT6l) . Then 

Var(Z^iv) = -didid2d2QN ■ 

Theorem 12.41 says that 
y8iT{ZgN,ip) = { — didid2d2Q N , (pM(p) = / QN{z,w)idd(p{z) Aidd(f{w) , (23) 

J MxM 

for all real [m — l,m — l)-forms on M with coefficients. 
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Since Pjj G C°°(M x M) and Pn{z, w) < 1 for z ^ w, for sufficiently large (so that the 
Kodaira map for is an embedding), it follows from (120!) that Qat is C°° off the diagonal, 
for > 0. 

Proposition 12.21 implies that the pluri-bipotential decays rapidly away from the diagonal: 
Lemma 2.5. |SZ4l Lemma 3.4] For b > y/j + q + 1, j > 0, we have 



\V^Qr,{z,w)\ = 0{^) , for d{z,w) > ^'^^ 



,A^V VN 
Proposition 12.31 yields the near-diagonal asymptotics: 

Lemma 2.6. |SZ4l Lemma 3.5] For b E we have 

Q^(zo,^o + ^) =G(e-5H^) + 0(iV-i/2+-) , for \v\<b^/I^. 



Recalling fl?T]) . we write. 



F{\) := G{e-') = / log(l - e-'n ds (A > 0) , (24) 



oo 



so that Qisi = F o {— logP/v). By Proposition 12. 3[ 

-logPjv (^^0,^0 + ^^ = ||^^r + 0(|t;pAr-i/2+^) < b^logN 



(25) 



It follows from Lemma [2.61 and fl25|) that Q^r G C^{M x M) and the first partial derivatives 
of Qjv vanish along the diagonal in M x M, for ^ 0. (We note that Qjv is C°° off the 
diagonal, but is not at all points on the diagonal in M x M, as the computations in jSZ4] 
show.) We furthermore have the near-diagonal asymptotics: 

Lemma 2.7. |SZ4[ Lemma 3.7] There exist a constant Cm G (depending only on the 
dimension m) and Nq = No{M) G Z"*" such that for N > No, we have: 

i) The coefficients of the current did2QN are locally hounded functions (given by point- 
wise differentiation of Qjsi), and we have the pointwise estimate 



\did2QN{z,w)\<CmN for 0<\w-z\<b^^-^. 

ii) Ifm>2, the coefficients of the current didid2d2QN ore locally U^~^ functions, and 
we have the estimate 



\didid2d2QN{z,w)\ < ^""^ for < \w - z\ < bx/^^^;^. 

\w — z\'^ \ N 

Lemma 2.8. |SZ4l Lemma 3.9] For N sufficiently large, 
-didid2d2QNizo,zo+^) = NVaTZiv)+0 {\v\~^ N^/^+') for < \v\ < 6^1^, (26) 
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where Var^ G T*,^'H^) ® V'^^^C^) is given by 

Var^(i;) := -^F^^\^\v\^) {v ■ dz){v dz){v ■ dv){v dv) 

- lF'^^\l\v\^) [{dz ■ dz){v ■ dv){v ■ dv) + {v ■ dz){v ■ dv){dz ■ dv) 

+ {v ■ dz){dz ■ dv){v ■ dv) + {v ■ dz){v ■ dz)(dv ■ dv)'j 

- \F"{^\v\^)[{dz-dv){dz-dv) + {dz-dz){dvdv)]. (27) 

Differentiating <^^, we note that 

^"(^) = A^A^' ^^'^(^) = -Acsch^A, F^^\X) = \ cothAcsch^A . (28) 



Thus, 



F(^)(A) = (-l)^^ ^-^—^ A-^+i + 0(1) (A > 0) 



for j > 2, and hence 



Var^o^^^_ J 0{\v\-') for \v\ > 

Var^(t;)-<; o(|^|4e-H^) for > 1 " ^^^^ 



3. The sharp variance estimate: Proof of Theorem 11.11 

3.1. The codimension one case. To illustrate the basic ideas of the argument, we begin 
with the proof for the case k = 1. By Theorem \2A\ we have 

Var(Z,^,^)= [ (z) iddifiiz) , (30) 

where 

I^{z)= [ QNiz,w)iddip{w) . (31) 

J{z}xM 

We let 

ml 

denote the volume form of M, and we write 

idd^ = i)VLM, ipeC^iM), (32) 

so that 

I''{z)= f QN{z,w)tlj{w)nM{w) . (33) 



J{z}xM 

To evaluate I^{zq) at a fixed point zq G M, we choose a normal coordinate chart centered at 
zo as in ^ and we make the change of variables w = zo + By Lemma [2751 and fl32|) - fl33|l . 

we can approximate T'^ {zq) by integrating fl33l) over a small ball about zq: 



rNi I ^ I . V \ ^ I V 



\v\<b^yogN 



1 {zq) = Qn { Zo, Zo + ^ ] ij { Zo + ^ ] { Zo + ^ ] + O 



(34) 



where b = + 2. 
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Since u: = ^ddloga = ^dd [\z\'^ + 0(|zp)] in normal coordinates, we have 

Sjk + o 



a; ^0 + 



2 

V 



ldv,Advk = ^dd\v\^ + 



, (35) 



for \v\ < by/log N. Hence 



N J ml 



(H _ 



\ogN 



(36) 



for |t>| < by/\ogN, where 



1 ^i 



^E{v) = ^i-dd\v\') =ll-dv,Adv, 



m m 



denotes the Euchdean volume form. Since (p E and hence '^(^ + :^) = i'iz) + 0{\v\/^/N), 
we then have by Lemma [2.61 and fl34p - (l5^ . 



AT" 



x{nE{v) + 0{N-'/^+')} 



G(e-2l-l') + 0(iV-i/2+e)| |^(^^) + 0(iV-i/2+.)| 
+ O (A^-™~^) 



V^(^o) 



lJ\v\<b^/]EgN 

Since G{e^^) = 0(e"^^) and hence 



G(e-5l-P)f]g(^;) + 0(Ar-^/2+^) 



G(e-5l-l^)r]^(t;) = OiN-""-') 



(37) 



(38) 



\v\>bylSiN 



we can replace the integral over the (6\/log iV)-ball with one over all of C™, and therefore 



■^{zq) 



(39) 



Recalling (I20]) . we have 



G'(e-^H^)(]^(^) 



47r2 



fc=i 



A;2 



1 °° r 
1 V ^ TT 



E 



TT 



m-2 



47r2 ^ fc^+s 

k=l 



Therefore, by ^ and (I39|)-(|l0l), 



Yax{Z,N,ip) 



M 



m-2 



C(m + 2) + 0(Ar-i/2+-'; 



C(m + 2) . 



^(;2)2fijvf(^) 



(40) 



(41) 



The variance formula of Theorem 11.11 for the case k = 1 follows from fl32|) and (jH 
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3.2. An explicit formula for the variance. In this section, we give an integral formula 
for the variance of simultaneous zero currents in higher codimension (Corollary 13.31) . which 
we shall use in the next section to derive the asymptotics of Theorem II. 1[ This integral 
formula is a modification of the following formula for the variance: 

Theorem 3.1. |SZ4t Theorem 3.13] Let 1 < k < m. Then for N sufficiently large, 
Var(Z,iv ... = did2 

where the current inside the brackets is an L} current on M x M given by pointwise multi- 
plication, is given by fl22l) . and EZgN is given by f|T3|) . Furthermore, \av(^ZgN^^ ,,n^ is 
an L} current on M x M if k < m — 1. 

By an current on the compact manifold M x M, we mean a current whose local 
coefficients are functions. 

Corollary 3.2. Let 1 < k <m. Then for N sufficiently large, 
Var(Z^]v^ gAf) = 81818282 

where the current inside the brackets is an current on M x M given by pointwise multi- 
plication. 

Proof. Let 

^ = ^[ ]Qn {-8idi82d2QNY'^ A {EZ,N K EZ.n)'"' 
j=i \3/ 

denote the expression inside the brackets in Corollary 13.21 regarded as a (4A; — 4)-form 
on M X M \ A, where A = {{z,z) : z G M} denotes the diagonal. By Lemma l2.7( ii). 
T = O (d{z, w)~'^^^^^ and hence T defines an current on M x M. It suffices to show that 
did2T is an current and thus 8182 can be moved outside the brackets in Theorem 13. 1[ 

Let = {{z,w) E M X M : d{z,w) < e} denote the e-neighborhood of the diagonal A. 
For a test form (f with m — k + 1 dzj^s, dzj^s, and dwj^s , and m — k dwj's, we have 

{d2T, ip) = - Mm T A 82^ = lim / c^aT A + lim / T A . 

Since T A ip = o(£:^^'"^^) on 8Us, the boundary integral goes to 0. By Lemma [2.7( ii) and 
the fact that Qn ^ C^{M x M), the pointwise-defined form c^aT is also 0{d{z, w)~'^^~^'^) and 
thus d2T is an current on M x M given by pointwise differentiation. 

Repeating the same argument with T replaced by and using part (i) of Lemma 12.71 as 
well as part (ii), we then conclude that did2T is an current. □ 



Y,(-iy ( .\di92QN A {8idi82d2QNT^ A {EZ,N m EZ.iv)'"'' 
.7=1 ^"'^ 



k 

.7=1 



Qn {-didid2d2QNY ' A (EZ,iv M EZ^n)" 



k-j 



Corollary 13.21 can also be shown directly, using the argument in the proof of Theorem 13.11 
in [SZ4j . 
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Corollary 3.3. The variance in Theorem \l.l\ is given by: 



E 



Qn {-didid2d2QNY ^ A {EZ,N m EZ.nY ^ A {dd^ K dd^) 



k-j 



where the integrands are in L^{M x M). 

3.3. Higher codimensions. Recalling ([HI), we write the formula of Corollary l3.3l as follows: 

k 

(42) 



Var([Z.....,.],^) = j:Qvfi^) 



A 



MxM 



Qn{z, w) [ - didid2d2QN{z, w)] 



LuizY-^ + O 
2k-2j 









m 


A 


+ ^ (^) 







M 



If A 



u^-^ A idd'^ + O 



where 



A idd(p{z) A idd(p{w) 
(43) 

1 



uiwY-^ Aiddif{w) + { — 



QNiz,w) [- didid2d2QNiz,w)\-' 'A 

'{z}x]\I 

e T*J-^^^-^{M) . (44) 

The integrand in fH4|) is regarded as an (m, m)-form (in the u; variable) with values in 
y*i-ij-i(]Vf). 

Fix a point ^ ^ ■, and let 1 < j < A; < m. To evaluate If [zq)^ we write 



cj'^-^' A I 



dt;-^ A dv^ 



\J\ = \K\ = m — j + 1 



By Lemma 12. 5^ we can replace integration over M in f l44l) with integration over the small 

i-i 



ball of radius 6 A/log N/N, with 6 = i/m + 4, to obtain: 



|?)|<6v^Iog7V 



<5Af Zq, Zo + 



N 



-didid2d2QN Zo, Zo + 



dv^ A dv^ + O {N-"^-^) . 



By Lemma [2.71 the above integrand is L^, and hence by Lemma [22 

i-i 



j^2j-2-m 



Fil\v?){Ys.v'liv 



y^iljjK{zo)dv^ A dv^ + 0{N^^ 

|t;|<6v'Iog7V ■' 

V^ji^(^o) / F{l\v\'')\yar'^{v)y~'dv' Adv"" + 0iN--2+n 



.(45) 
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Here, we replaced the integral over the (6v^IoglV)-ball with one over all of C™, since by fl2I?l) 
we have F(||f p){Var^(t')}-'^^ = 0(e~l''l^) for |f | > 1, and hence 

1 J-i 



/ 

J\v\ 



'\v\>by/YS^ 

It follows from (USD and (USD that 



Jm 



(46) 



where = {Bj^^^} is a universal Hermitian form on T*"^ J+i(M). Theorem 11.11 

then follows from (1421) and (1461) with 

k 

□ 

4. Asymptotic normality: Proof of Theorem 11.21 

The proof of Theorem 11.21 is an application of Propositions I2.2H2.3I to a general result of 
Sodin-Tsirelson [ST] on asymptotic normality of nonlinear functionals of Gaussian processes. 
Following [STJ, we define a normalized complex Gaussian process to be a complex- valued 
random function w{t) on a measure space {T,fi) of the form 



w 



where the Cj are i.i.d. complex Gaussian random variables (of mean 0, variance 1), and the 
Qj are (fixed) complex-valued measurable functions such that 

l^i(^)r = 1 for all t e T. 

We let w^^w'^^w^, ... be a sequence of normalized complex Gaussian processes on a finite 
measure space (T,ix). Let /(r) G L'^(W^,e~^ ^'^rdr) and let ■?/' : T ^ M be bounded measur- 
able. We write 

^^(^^)= [ /(k^(t)|)^(t)rf/i(t). 



Theorem 4.1. [STl Theorem 2.2] Let pAr(s,t) be the covariance functions for the Gaussian 
processes w^{t). Suppose that 

• ^ .^^ J^J^\pN{s,t)\^''ij{s)ij{t)d^{s)d^{t) ^ ^ 

N^cc sup, J^\pNis,t)\d^M{t) 

for a = 1 if f is monotonically increasing, or for all a G Z+ otherwise; 
ii) lim sup / \pNis,t)\dfi{t) = 0. 

N^OO gfzT Jrp 

Then the distributions of the random variables 



Var(4) 
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converge weakly to A/'(0, 1) as N ^ oo. 

We apply this result with /(r) = logr and (T, /x) = {M,Qm)- To define our normalized 
Gaussian processes on M, we choose a measurable section ai '■ M ^ L oi L with 
= 1 for all z G M, and we let 

be an orthonormal basis for H^[M, L^) with respect to its Hermitian Gaussian measure, for 
each iV G Z+. We then let 



9j {z) := — = ; , j = l,...,d7v. 



Since = ||S'j^||/jjv, it follows that = Yli^jdf defines a normalized complex Gauss- 
ian process, for each N G (where the Cj are i.i.d. standard complex Gaussian random 
variables). In fact, 

where 

is a random holomorphic section in H^{M, L^). The covariance functions Pn{z, w) for these 
Gaussian processes satisfy 

\Pn{z,w)\ = Pn{z,w) . 
We now let (f he a fixed real (m — 1, m — l)-form on M and we write 



1 _ 

—ddif = ip VLm ■ 

71 



Then ip & C^, and 



M 



Zti^"") = I (log y'^mh- - log V^n{z,z)^ ^ddifiz) = {Z,N, ^) + 



where the k^ are constants (depending on L — > M and ip, but independent of the random 
sections s^). Hence Zff{w'^) has the same variance as the smooth linear statistic (Z^jv,*/?), 
and it suffices by Theorem 14.11 to check that the covariance function satisfies conditions 
(i)-(ii) of the theorem. We start with (ii): by Proposition 12.21 

lim sup / PN{z,w)Q,Miw)=0- 

On the other hand, since Pn{z, w) < 1, it is obvious that the same limit holds for d{z, w) < 
verifying (ii). 



To check (i), we again break up the integral into the near diagonal d{z, w) < b\/ and 



the off-diagonal d{z, w) > b^^^j^. As before, the integrals over the off-diagonal set tend to 
zero rapidly and can be ignored in both the numerator and denominator. 
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On the near diagonal, we replace Pn by its asymptotics in Proposition [231 The asymptotic 
formula for Pn has a universal leading term independent of z and has uniform remainder, 
so condition (i) (with a = 1) becomes 

limmf — — J > 0. 

Since ip E C^, the ratio clearly tends to 2"'" J^^ "ipi^Y Qm > 0, which verifies (i) and completes 
the proof of Theorem ll.2[ □ 

5. Open problems on smooth and counting statistics 

In this section, we present a number of open problems on smooth and discontinuous linear 
statistics of zeros. 

(1) Asymptotic normality of the smooth linear statistics {Z^n s^j'P) only been 
proved in codimension one, i.e. when k = 1. But these random variables are likely 
to be asymptotically normal for all dimensions m and codimensions k. It would be 
interesting to prove (or disprove) this statement, in particular for k = m. 

(2) To our knowledge, no results are known to date regarding the asymptotic normality 
of the counting statistics A/'jv(s{', . . . , s^), even when m = 1. This is analogous to, 
but presumably harder than, the smooth linear statistic when k = m. 

(3) In |Zhj . Qi Zhong obtained surprising results on the expected value of the 'energy' 
random variable 

£g{s^) = '^G{ai,aj), Z^n = {ai, . . . , a^} , (47) 

summing the values of the Green's function G over pairs of distinct zeros of a random 
polynomial or section G H^{M, L^) of the A^-th power of a positive line bundle 
(L, h) over a compact Riemann surface M. (The Green's function is normalized to 
equal +oo on the diagonal.) Zhong proved that when G is the Green's function for 
the Riemannian metric induced by the curvature of /i, the expected energy has the 
asymptotics E£^G'(PAf) ~ -^^log^- It is known (N. Elkies) that -^A^logA^ is 
also the asymptotic minimum for the energy sum ( 1471) . The energy is a partially 
smooth linear statistic 

(G,Z,iv®Z,iv-Az^) (48) 

on M X M, where Az ^ represents the diagonal terms of Z^n ® Z^n. The statistic 
(HHl) is not smooth since G has a logarithmic singularity and since we subtracted the 
diagonal current. It is a random variable of one section sat in dimension one, but G 
is a function on M x M, so the variance of Sq involves the rather complicated four- 
point correlation function of Zg rather than the pair correlation, for which Theorem 
12.41 gives a useful formula. It would be interesting to investigate the variance of the 
energy Sq- It seems that it should tend to zero with N since the Gaussian measure 
is concentrated on 'polynomials' whose zeros are asymptotic minimizers. 

(4) The expected distribution of zeros can have quite disparate asymptotics when the 
ensembles are given Gaussian measures induced from inner products on the space 
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of polynomials (or sections) which use non-smooth volume forms or non-positively 
curved line bundle metrics. For instance, in the case where the measure is supported 
on an analytic plane domain C C or on its boundary, it was shown in [ SZ3j that the 
the expected distributions of random zeros of random polynomials of degree tend 
to the equilibrium measure of Q. This result was generalized to higher dimensions 
and more general metrics and measures in a sequence of papers |Bllt IBSt IBe] in 
which it is shown that the expected distribution of zeros tends to an equilibrium 
measure adapted to the measure and metric. In [Sh\ . an upper bound was given for 
the variances of the smooth linear statistics when the inner products are defined by 
arbitrary measures (and also for more general sequences of ensembles of increasing 
degrees). This bound is sufficient to prove that sequences of random zeros in these 
ensembles almost surely converge to their equilibrium measure, although the bound 
is not always sharp. The Szego kernels for the inner products in |SZ3t IBllt IBSj IBe] 
are quite different from Szego kernels for positive line bundles in this article, and so 
the asymptotics of the variances might be quite different. It would be interesting to 
determine them. 

(5) Results on expected values for the analogous SO (2) and SO(m + 1) ensembles of 
random real polynomials in one or several variables were given by [BD[ IEK[ IRoj ISSj , 
and a (global) variance result for real zeros was given in [WsJ. Maslova [Ma] proved 
the asymptotic normality of the number of real zeros for the Kac ensemble |Kac] of 
random real polynomials on (as well as for some non-Gaussian ensembles). But as 
far as we are aware, asymptotic normality for numbers of real zeros in the SO(m + 1) 
ensemble has not been investigated. 
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